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ABSTRACT 

Radiative deceleration of relativistic jets in active galactic nuclei as the result of 
inverse Compton scattering of soft photons from accretion discs is discussed. The 
Klein-Nishina (KN) cross section is used in the calculation of the radiation force due 
to inverse Compton scattering. Our result shows that deceleration due to scattering 
in the KN regime is important only for jets starting with a bulk Lorentz factor larger 
than 10"^. When the bulk Lorentz factor satisfies this condition, particles scattering 
in the Thomson regime contribute positively to the radiation force (acceleration), 
but those particles scattering in the KN regime are dominant and the overall effect 
is deceleration. In the KN limit, the drag due to Compton scattering, though less 
severe than in the Thomson limit, strongly constrains the bulk Lorentz factor. Most 
of the power from the deceleration goes into radiation and hence the ability of the 
jet to transport significant power (in particle kinetic energy) out of the subparsec 
region is severely limited. The deceleration efficiency decreases significantly if the jet 
contains protons and the proton to electron number density ratio satisfies the condition 
np/neo > 27inin/Mp where 7min is the minimum Lorentz factor of relativistic electrons 
(or positrons) in the jet frame and fj,p is the proton to electron mass ratio. 
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1 INTRODUCTION 

High energy observations of blazars strongly suggest that 
, the emission comes from relativistic jets of active galactic 
' nuclei (AGN) (e.g. von Montigny et al. 1995; Thompson et 
' al. 1995; Gaidos et al. 1996; Quinn et al. 1996; Schubnell et 
. al. 1996). Other evidence for relativistic AGN jets includes 
' the observation of variability in radio, optical, and high en- 
ergy gamma ray emission from blazars. Although there is 
no consensus on the details of AGN, it is widely accepted 
that a jet may form near a black hole with an accretion disc 
(e.g. a review by Blandford 1990). Since a substantial frac- 
tion of the binding energy of accreting material is dissipated 
and converted to radiation, the disc is a strong source of 
soft photons (e.g. as black body radiation or reprocessed ra- 
diation). Interaction of relativistic particles in the jet with 
these surrounding photon fields can be important, and may 
contribute to the observed gamma ray emission, e.g. through 
inverse Compton scattering (Reynolds 1982; Melia & Konigl 
1989; Dermer & Schlickeiser 1993). Disk emission can also 
affect the jet dynamics near the black hole through radiative 
acceleration or deceleration. The effect of Compton scat- 
tering on the jet flow was discussed by O'Dell (1981), and 
Phinney (1982, 1987), and it was suggested that under cer- 
tain (quite stringent) conditions, the jet can be accelerated 
through Compton scattering. However, as pointed out by 
several authors (e.g. Phinney 1987; Melia & Konigl 1989; 
Sikora et al. 1996), in the region close to the black hole. 
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Compton scattering can be more effective in slowing down 
rather than accelerating the jet, and can constrain the bulk 
flow of the jet in that region. 

There are AGN emission models in which protons are 
accelerated to ultra high energies (e.g. Konigl 1994). One 
of the attractive features of this type of model is that pro- 
tons can be accelerated to high energies without significant 
energy loss. These high energy protons can interact with 
photon fields producing electron-positron pairs. It is possi- 
ble that through pair cascades, an ultrarelativistic jet is pro- 
duced. Alternatively, an ultrarelativistic jet can be produced 
through rapid acceleration of such as by magnetic re- 
connection (e.g. Haswell et al. 1992) or by rotation-induced 
electrostatic potential drop as in pulsars (Michel 1987). In 
both cases, the bulk flow of the jet may initially have a large 
Lorentz factor, and plasmas in the jet can be highly relativis- 
tic. Therefore, scattering in the Klein-Nishina (KN) regime 
is important, and should be included in calculation of radia- 
tion force due to inverse Compton scattering. There is then 
the question of whether the jet with initially large bulk flow 
speed is subject to the same (severe) radiation deceleration 
as in the Thomson regime, and this will be explored in this 
paper. 

The constraint on the bulk flow by Compton scattering 
should strongly depend on the soft photon distribution of 
the disc. In our discussion, disc emission is assumed to be 
axisymmetric, and so the only effect on the jet flow consid- 
ered is in the jet direction (it is usually assumed that the 
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jet is normal to the disc plane). However, there are cases in 
which accretion may not be circular, e.g. eccentric accretion 
can occur if the black hole is in a binary system (e.g. Begel- 

rnan, Blandford & Rocs 1980; Eraclcous ct al. 1995), and 
this possibility will be considered in a forthcoming paper 
(Luo 1998). 

The deceleration efficiency strongly depends on the 
composition of the jet. An AGN jet that contains protons 

would be subject to less severe deceleration because protons 
scatter photons with much smaller cross section than that 
for electrons (or positrons) and this effectively increases the 
jet inertia. 

In this paper, we extend calculations of the radiation 
force by including the Klcin-Nishina effect as the earlier cal- 
culations were done only in the Thomson scattering regime 
(e.g. O'Dell 1981; Phinney 1982; Sikora et al. 1996). In Sec. 
2, the average radiation force on the bulk plasma flow is de- 
rived in both the Thomson and Klein-Nishina limits. Con- 
straints by radiative deceleration on the Lorentz factor of 
the jet are discussed in Sec. 3. 



2 RADIATIVE EFFECT ON JET FLOW 

Through Compton scattering of external photons, individual 
particles in plasmas lose energy and at the same time there 
is momemtum transfer to the plasma. For plasmas in a jet, 
the momemtum transfer can affect the jet dynamics, i.e. 
the bulk flow can be either accelerated or decelerated (i.e. 
radiative drag). 

Consider a comoving cell with energy E (in nieC^ per 
unit volume) in a relativistic jet with the bulk Lorentz fac- 
tor r = 1/(1-/3^)^^'^ where is the bulk velocity in c. 
The tilded quantities correspond to the values in the lab 
frame (K). We assume that the jet contains mainly electron- 
positron pairs (there may be protons as well but in small 
number) and that within the cell electron-positron pair pro- 
duction is not important. The latter assumption may not 
be accurate since pairs can be injected continuously along 
the jet. Nonetheless, it allows us to single out the radiative 
effect due to Compton scattering, and the result should not 
change qualitatively when the pair production effect is in- 
cluded. Let / be the radiative force (in meC per second per 
unit volume) on the coll seen in the jet comoving frame (Kj). 
Then, the rate of change of T due to radiation force on the 
cell is given by (e.g. Phinney 1982; Sikora et al. 1996) 

/3b (2.1) 



E 



where E is the cell energy in Kj, and we assume that the jet 
is along the z direction. Then, for fz > 0, the jet is acceler- 
ated away from the source, and for fz < 0, it is decelerated 
toward the source, i.e. radiation drag. In (2.1), the average 
force fz depends on the particle distribution which is influ- 
enced by the specific acceleration or energy loss mechanisms. 
In the following discussion, we only consider radiative effects 
due to Compton scattering. 



the accretion disc emission, and is axisymmetric with respect 
to the jet axis (which is in the z direction). A schematic 
diagram for such a disc-jet system is shown in Figure 1. 

Emission from the disc is modeled as the sum of emission 
from series of rings centered at the black hole; each of them 
emits photons with a characteristic energy e^^ where 7? is the 
radius of the ring. The flux number density of photons from 
the ring Rr^ R + dR is given by F(e, R) = F{R)S{£ - e^) 
where e — r(l + /?b cos6)e is the photon energy (in meC^) 
seen in K, e is the corresponding energy seen in Kj, 9 is the 
angle between the photon propagation direction and the jet 
axis, and F{R) is given by (Blandford 1990) 
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where 7 = (1 — Rmin/R)^^^, Rmin is the radius of the inner- 
most part of the disc, M is the mass of the black hole, G 
is Newton's constant, I/Edd is the Eddington luminosity. We 
assume i?max to be the maximum radius of the disc, within 
which emission can be adequately described by (2.2). 

A detailed model for emission from the accretion disc 
should include effects such as reprocessing of radiation, disc 
corona, etc. As we are concerned with the radiative effect on 
the jet flow, all these details will be ignored and we assume 
that disc emission is blackbody with e^j w 2.7k^T[R) /me(? , 
where T{R) is the effective temperature of emission from 
the ring with R. Then, for an optically thick disc, T(R) 
is identified as the surface temperature, and it is given by 
T{R) = {SGMMI/STTR^assy^*, that is. 
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where R > Rmin, tsb is the Stefan-Boltzmann constant. 
Figure 2 shows plots of F{R) for Rmin/Rg = 5 and 6, where 
Rg = GM/c^ « 1.48 X lO"cm(M/lO*M0) is the gravita- 
tional radius. The fiuxes are peaked at Rq — {7/G)R^i^ ~ 
5.8Rg and 7Rg, respectively. 

Let nph(e, f2, i?)rf7? be the number density of photons 
emitted from the ring between R and R + dR, where s and 
fl — {<j),cos6) are the energy and direction of incoming pho- 
tons, respectively. Using (2.2), the photon number density 
(per unit R), nph(£, Q, R), in Kj can be written as (e.g. Der- 
mer & Schlickeiser 1993) 



nph(e,f2, i?) = 



RF{£, R) 
27rc(i?2 + z2 
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where cosfii^ = z/{R^ -\- z^Y^"^ , z is the distance from the 
disc surface, is the angle of incoming photons relative to 
the the jet axis in Kj . 



2.1 Radiation from the disc 

To study the radiative effect on jet dynamics due to inverse 
Compton scattering, a specific model for the external photon 
field is required. We assume that the external photon field is 



2.2 Radiation force 

The radiation force / (on the jet bulk flow) can be derived 
by calculating the rate of average momemtum transfer to 
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Figure 1. A schematic diagram of an AGN jet. 
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Figure 2. Photon flux density from the ring as a function of 
R. We assume M = IO^Mq, M = O.lLEdd/c^- The upper and 
lower curves are obtained by assuming iiniin = 5ilg and GRg, 
respectively. 

plasmas in the jet as a result of inverse Compton scattering, 
which is given by 

f = / dQed'y ^ ne(7, fie) 



(2.5) 
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where Z) = l- /3-k, 7 = 1/(1 - (i^y^"^ is the particle 
Lorentz factor in the jet frame (Kj), n^^^{e, Q.) is the photon 
density, p = 7/3 is the momemtum (in mec), ne(7,ne) is 
the particle distribution, Qe ~ (<^e,cosSe) is the direction 
of electron motion, da/desdO,s is the differential cross in 
Kj whose approximations in the Thomson and KN regimes 
are given respectively by (Al) and (A2). The direction of 
incoming and scattered photons in Kj are represented by k 
and ks, respectively. All quantities with subscript s are for 



scattered photons. As we use axisymmetric disk emission as 
given by (2.2), the only relevant component of the force is 
along z axis. 

We first derive the radiation force in the Thomson limit. 
In the Thomson scattering regime, scattering is elastic and 
we have e'3 w e' where e' and are respectively the energies 
of the incoming and scattered photon in the electron rest 
frame {K'). Further, we assume a beaming approximation 
in which the scattered photons propagate approximately in 
the electron's direction in Kj provided that 7 1. Based 
on these considerations, an approximation for the differential 
cross section in Kj can be obtained as (Al), and the force 
is calculated as 



=o-T / 
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where we assume that the disc emission is important within 
Rmin < R < Rmux, cT^ Is the Thomsou cross section, the 
average is made over the electron (or positron) distribution, 
(...) = J ciOed7 (...)ne(7, f2e)/n-eO, and rieo is the average 
number density. The distribution ne(7, fie) is assumed to be 
isotropic in Kj with a power law: 



ne(7,0e) = ^"e07 



(2.7) 



where p is the electron spectral index with values between 



2 and 3, 7min < J < Jir 



the constant Co is chosen 



such that (2.7) is normalized to rieo- We assume the cy- 
clotron time scale is much less than the relevant time scales 
of acceleration or deceleration, electron (positron) energy 
loss. Therefore, all quantities in (2.7) can be regarded as 
average over the cyclotron time. For 7^^,,^ 3> 7„i„ and 
p > 2, we have Co ~ (p — l)7m7ii^- Using (2.7), we have 
(7^) ~ [{P - l)/(3 - P)]7^7„'(7i;f- 7ii7i') which reduces to 
(7') « 272 ^„ ln(7.„.x/7.„i Jfor7= 3. °"° 

It follows from (2.6) that acceleration (/^ > 0) occurs 
for cos > /3^^ corresponding to when electrons see most 
incoming photons with 6*, ~ 0, and deceleration (/^ < 0) 
for cos^f^ < /9j, corresponding to when electrons see most 
photons with 9r m tt in the jet frame. Therefore, the key 
role in radiative acceleration or deceleration is played by 
anisotropic scattering, which can be understood as follows. 
Assume a cell moving along the jet containing relativistic 
plcisma with an isotropic distribution. The particles in the 
jet frame see anisotropic photon fields, i.e. the radiation 
mainly comes from behind or in front depending on F of the 
bulk motion. When cos 0, < P^, in the jet frame, although 
each electron would scatter photons into its direction of mo- 
tion, on average there is more momemtum beamed along 
the jet direction. Thus, the cell is subject to a force in the 
opposite direction. The critical F, above which the radiative 
deceleration (/^ < 0) occurs, can be derived from (2.6) by 
setting /z = (e.g. Sikora et al. 1996). Since the critical F 
is generally small (e.g. Phiney 1987; Sikora et al. 1996), our 
discussion concentrates on deceleration of jets with a large 
initial F. 

In the case of deceleration, since the flux (2.2) is peaked 
aX R = Ro = (9/7)i?„i„, an approximation for fz at z ^ Ro 
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is given by 



RgC / Ro \ * 2 



(2.8) 



where /Up = mp/me is the proton to electron mass ratio, 
Ro = 7R^ia/6 is the radius at which the flux is peaked, 
7] = Ld/L^^^ is the radiation ofRcicncy of disc emission, 
= 2'KmeC^ J dRRF{R)£i^ is the luminosity of the source, 
and the equality is for p = 3. The radiation force increases 
rapidly with . 

Although (2.6) is derived using the beaming approxima- 
tion as described by 5{Qs — (^e) in Eq. (Al), apart from /3 ~ 1 
in the average (7^/?), it is in good agreement with the result 
derived by Sikora et al. (1996) using covaxiant formalism. 
One may also reproduce the result given by O'Dell (1981) 
where the photons are from a point source, for which we may 
assume 9^=0 and that 27r / dRRF(R)e„ L^JnieC^ is 
the luminosity of the source. In this special case, the parti- 
cles are accelerated away from the source. 

In calculating (2.6) and (2.8), we assume e' < 1. In 
the electron rest frame, the photon energy is e' = e'^ = 
7(1 — /3cos©H)r(l — cos9j^)eR, where Or is the angle 
between the incoming photon direction and electron motion 
(in Kj). Then, the condition for all particles to scatter in 
the Thomson regime is 



7e« < 1/2. 



(2.9) 



with — r(l — /9(, cos 6'_f?)ej, . This condition requires that 
photons of the highest energy, seen by particles with = 
TT, satisfy e'^ < 1, which may not be the case if F, or the 
average Lorentz factor ((7)) of the plasma in Kj, are large. 

2.3 The Klein-Nishina regime 

For an ultrarelativistic jet with a large bulk Lorentz factor 

or when particles are highly relativistic in the jet comoving 
frame, Compton scattering should be treated in the Klein- 
Nishina scattering regime. The minimum F such that all 
particles of energy 7 are in the KN regime can be derived 
as follows. In the jet frame, since the soft photon direction 
is 6r « TT, electrons with 9e = n see the lowest energy 
photons. Since cos Qr — cos 6r cos 6e + siuOr sin 6e cos(</) — 
06 ) ^ — cos^^cos^e (where ([> and 0e are the azimuthal 
angles of the photon and electron), and from e'^ > 1, we can 
derive the condition for all particles to scatter in the KN 
regime. 



> 



(2.10) 



(2.11) 



7 eR{l- P^cosOr) 
or 
F 

- < +5080^), 

where we assume that the plasma is relativistic in the jet 
frame. For er = 6 x 10"'\ we have F/7 > 5 x 10''/(1 - 
cos0r) or F/7 < 3 X 10"'^(1 -I- cos^fl). Inequalities (2.10) 
and (2.11) are the KN condition for particles of 7 with any 
9e , which requires an extreme ratio (very large or very small) 
F/7. However, in the following discussion we show that the 
condition for the KN scattering to be the dominant process 
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Figure 3. Plots fic as a function of R/Rg for different z and 7. 
The solid and dashed curves correspond respectively to 7 = 10^ 
and 10^. In each case, the curves from bottom to top correspond 
to z/Rg = 5, 10, 30, 50. We assume the bulk Lorentz factor 
r = 10*. 
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Figure 4. The minimum T for KN scattering to be dominant in 
deceleration are plotted as functions of R/Rg. The curves from 
top to bottom correspond to z/Rg = 20, 10, 5, 2. We assume 
Tmin ~ ^min = 5Rg, and use the temperature function (2.3) 
with M = O.lLEdd/c^- 

in radiative drag is less strict than (2.10), (2.11), and is 
generally satisfled for moderate F. 

In general, given a particle distribution, the Thom- 
son and KN approximations apply respectively to parti- 
cles with cos^e < fic and cos6'e > /itc, where Hc satisfies 
F7(l — P)j,c cos Or){1 — Pi cos 9R)eR = 1, that is. 



- r7(l - Pb COS0r)£r 

TyiPb - cosOR)eR 



for 7£„ > 1/2; 
for 7£h < 1/2. 



(2.12) 



Thus, the conditions (2.10), (2.11) correspond to an extreme 
case; if one of these two conditions is satisfied, then fic = 
— 1, and in this case the KN approximation applies to all 
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particles. Figure 3 shows plots of /Xc as functions of R/Rg. 
For given z and 7, the KN scattering regime is located above 
the curve (i.e. for electrons or positrons with cos^e > A*c). 

To calculate fz, particles of energy 7 in the cell are di- 
vided into two parts according to (2.12). Let /z,t and /z,kn 
be the force due to scattering by particles with cos Be < A*c 
(the Thomson regime) and cos^e > Hc (the KN regime), 
respectively. Then, the radiation force can be expressed as 



fz = fz,T + /z,KN. 



(2.13) 



where fz,T and /z,kn can be calculated approximately using 
the cross sections (Al), (A2) as given in the Appendix (cf. 
Eq. A4 and A5). Assuming that incident photons in the 
jot frame arc approximately beamed, cos^^ ~ — 1, which is 
valid for r S> 1, we have 
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1 
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(2.15) 



where Hi-^e^ - 1/2) = 1 for 72^ > 1/2 and H{-fe„ - 1/2) = 
for 7e^ < 1/2. In the Thomson approximation, /ic = 1, 
we have fz = fz,T, which reduces to (2.6) with cos 0^ w — 1. 
For He < 1, fz,T can be positive even though we assume 
cosOjf = —1. This is because /z.t includes particles with 
— 1 < cos^e < He and if most of them have 9e > 7r/2 we 
would have /z_t > 0. In contrast, /z,kn is always negative 
since it includes only those particles with fj,c < cos9e < 1, 
most of which move forward along the jet and contribute to 
drag. From (2.14) and (2.15), we sec that when /ic < 1 scat- 
tering in the KN regime becomes important while the force 
due to scattering in the Thomson regime starts to decrease. 
From (2.12), this condition becomes 

1 



r> = . (2.16) 

27(1 -/3b cos6l^)e„ 

Figure 4 shows the condition (2.16) as a function of R with 
7 being replaced by (7). We assume 7 = 80, p = 3. Thus, 
the average Lorentz factor in the jet frame is (7) = 160. 
This gives, for example, P > 400 for z = 2Rg, and T > 700 
for z = 5Rg. These lower limits can be further reduced as 
scattering in the KN regime extends to higher 7. 

Integration over 7 in Eq. (2.14) and (2.15) can be done 
analytically using (2.7), and is given by Eq. (A8), (A9) in 
the Appendix. For 2ejj < l/7„ij,, (2.15) has the following 
approximation 

/z,KN « 



— aa^Ueo I 



dR 



RF{R) 3(£^7„.J 



(2.17) 



-/z,KN/<7>°eO 
-.^,T/<7>°eO 




Figure 5. Variation of /z,T/(7>"e() (dashed), /z,KN/(7>"'eO 
(solid) with r for z = 5Rg (upper) and z = lORg (lower). The 
dotted curve (for z = 5Rg) is calculated by ignoring scattering in 
the KN regime. 

with a — 2/(p— 1) + ln(64e)/p — In 2/(p + 1), which ranges 
from 4.4 (p = 2) to 2.6 (p = 3). For large F satisfying (2.16), 
fz,KN overtakes fz,T and become important. For even larger 
F, we find 



fz^fz 



1 -I ^ 



dR 



RF 



7„i„P 



27r(i?2 + «2) 
ln(2£„7„in) 



(2.18) 



where Ri is the radius of the ring which emits photons with 

energy e„ = 1/(27^;J. If > l/(27,„iJ for < R< 

^maxi ^1 is replaced by -R„i„. Figure 5 shows fz,T / {'y)neo 
(Eq. 14), fz, KN / {'Y)neo (Eq. 15) as functions of F for z = 
5Rg, lORg. The disc parameters are the same as in Figure 
2. We assume p = 2.8, = 5Rg, = 507?g. The 

value of fz is significantly lower than that calculated us- 
ing the Thomson approximation (dotted curve), indicating 
the importance of the KN scattering. The KN scattering 
starts to overtake the Thomson scattering when F > 300 for 
z = 5Rg and F > 200 for z = WRg, and it is described by 
Eq. (2.17). Scattering in the KN regime becomes dominant 
in decelerating the jet at F > 10^ and the approximation 
(2.18) applies. The force component due to scattering in the 
Thomson regime becomes positive for large F since most 
particles that scatter in the Thomson regime move toward 
the black hole and hence contribute to acceleration of the 
cell. However, the net force is still negative because the com- 
ponent due to forward moving particles that scatter in the 
KN regime dominates. For extremely large F, we may have 
2£^ > P7„i„/(p— 1) lii(2e^7^;j,), and then, the second term 
on the right hand side of (2.18) becomes important, and 
fz *i /z,KN increases logarithmically with F. 



3 RADIATIVE DECELERATION 

A relativistic jet can form in the vicinity of a black hole 
through cither (1) a hydrodynamic process, in which ra- 
diation from an accretion disc plays an important role in 
producing and coUimating the jet, or (2) a MHD process. 



6 Luo and Protheroe 



in which magnetic fields play major role in the jet flow. In 
either case, regardless of the details of jet production, which 
are poorly understood, the jet must pass through intense ra- 
diation fields from disc omission and be subject to radiative 
deceleration. In the following discussion, we assume that the 
jet has initially a large F and then undergoes deceleration 
due to the radiation force discussed in earlier sections, and 
we examine the constraint on the bulk flow when scattering 
in the KN regime is included. 



3.1 Constraint on F 

The bulk Lorcntz factor F as a function of z can be cal- 
culated by integrating (2.1) along the jet direction. For de- 
celeration from the KN to Thomson regime, we use (2.13) 
together with (2.14) and (2.15) for which integration over 
R can be done numerically. Wc first consider an electron 
positron jet with an energy density (in niec) given by 
E = rieoi'y}- (A relativistic jet containing protons will be 
considered in Sec.3.2.) Prom the distribution (2.7), we have 



'p-2 



p-2- 



7„ 



7„ 



(3.1) 



for p > 2 and 7max 3> 7min. Figure 6 shows dlnV/dz 
as a function of z and 7. Radiative drag is severe, i.e. 
dhiT/dz <C —1, for small 2 and moderate F, but decreases 
as z increases. Plots of dlnF/dz as a function of z for an 
electron-positron jet are given in the solid curves in Figure 
7. For a very largo bulk Lorontz factor, the deceleration ef- 
ficiency decreases significantly because scattering is in the 
KN regime. For example, we have a much smaller d\nT/dz 
for F = 10® , as shown in Figure 7. 

In the Thomson limit, we can solve (2.1) with (2.8) to 

find 



r{z) 



(3.2) 



(3.3) 



1 + V^O^T ' 

where Fo is the initial value at z = zo, and 
For p = 3 and the distribution given by (2.7), reduces to 

(r^^o'!t{^)'{^-i)-r..J''i-r...h..l (3.4) 

Assuming that the initial value of F is large and 
r;FoCT > 1> 



one has the asymptotic value given by 
F(oo) ^ 



r)^^{zo,oo)' 



x/7„i„ = 5 X 10^ 7„i„ = 50, 
ig, even with Fo 3> 1, we have F(cx)) 



(3.5) 



2Ra 



1 for 



For, p = 3, 7„ 
and Zo = lORg 
rj = 0.1, and F(c») ^ 10 for rj = 0.01. Thus, in the Thomson 
limit, radiation drag is indeed a severe constraint on F as 
already shown in earlier work, e.g. by Melia & Kongl (1989), 
Sikora ot al. (1996). 

Deceleration described by (3.2) applies only for small F 
or 7^^^. When F satisfies (2.16), deceleration is mainly due 



log(-dr/rd^) 




Figure 6. A plot of -cflnP/ciz vs. z (with z = z/Rg) and P. The 
disc pajrameters are the same as in Figure 5. 




Figure 7. Plots of —dlnF/dz vs. z (z = z/Rg). The solid curves, 
from bottom to top corresponding to V = 10''', 50, lO'^, 10'^, are 
for a pm'c electron-positron jet (i.e. Up/n^a = 0, cf. Eq. 3.1). Tlie 
dotted and dashed curves are for tlie jet containing cold protons 
with rip/rieo = 0.1 and 0.5, rcsijcctivcly (cf. Sec. 3.2). In each 
of these two cases, the curves from bottom to top correspond to 
P = 50, 10^, 10^. The disc parameters axe the same as in Figure 
5. 

to particles scattering in the KN regime. From (2.17), and 
using (2.1), the bulk Lorentz factor is calculated to be 

Fo 



F(2) 



with 



2)- 



3q(p-2)X£T fls 
167r2P-i(2p- 1) 7?o 



2p-l 



(f) 



2p-l-| 



7 

' rr 



(3.6) 



(3.7) 



where p > 2, a = 4.4 — 2.6, e„ = 5 x 10~ . Deceleration 
is significantly slower than that calculated from the force in 
the Thomson approximation since is smaller than ^.p by 
a factor of 57„i„ lii(7max/7min) (fo'^ P = 3). Eq. (3.6) applies 



Relativistic jets 7 



only within a narrow range of F (i.e. 1/(7) < < l/7min) 
as shown in Figure 5. 

For even larger Fo, the drag is mainly due to scattering 
in the KN regime (2.18). Retaining only the first term on 
the right hand side of (2.18), we have 

for zo < z < zi, where zi is the distance at which r{zi) is 
so small that (3.8) docs not apply. Figure 8 shows the rapid 
drop of T{z) with the distance z, where wo assume that the 
jet starts with Fo = 10^ at zo = 5Rg and 10 Jig, respectively. 
In each case, the first curve from the left corresponds to 
the deceleration of an jot. The solid curves correspond 
to deceleration in the Thomson regime, and are calculated 
using (2.6). The bulk Lorentz factor decreases from F = 10^ 
to F = 10^ within a distance Az <^ zo. Radiation drag 
is effective except for Fo > (3?7^p/167r{7)e^)(i?g/2:o) ^ 5 x 
10'^{ri/0.1){Rg/zo) (for (7) = 100). However, for such a large 
Fo, the second term becomes important, and F decreases 
exponentially according to 

ln(F/Fo) « 

_ 3r?Mp(p- 1) Z"^, , i?„ 1 - cos6i 
Airp 

For p = 3, we have 

,„,r/r.,.-.3(i)i«(,- 1)1^1^, (3.0, 

with = Foe^. The drag strongly depends on 7„i„, and 
can be reduced for large 7„i„. 

3.2 A jet containing protons 

If a jet is heavily loaded with protons, the effect of radiation 
drag is reduced significantly, since protons scatter photons 
with much a smaller cross section compared to electrons, 
and the presence of protons effectively increases the inertia. 
When the jet contains cold protons, the energy density is 
E = 2neo7„i„ + "■pA'p where rip is the number density of 
protons and the electron spectral index is assumed to be 
p = 3 (cf. Eq. 2.7). Suppose that rip is so large that the 
condition 



(7)7mi 



(3.9) 



27min'^e0 



> 1 



(3.11) 



is satisfied. The value S in Eq. 3.11, corresponding to the 
ratio of the cold proton energy density to the electron energy 
density, can increase significantly if protons are relativistic 
in the jet frame. Then, fz/E is reduced by a factor 5. The 
asymptotic r(oo) given by (3.5) would increase by a factor 
S. For example, for 7^^^;^ = 50, one must include the inertia 
of protons if Hp/neo > 0.05, and the drag force would be 
significantly reduced. 

The effect of protons on the deceleration efficiency is 
shown in Figure 7 by dotted curves for Up/n^o = 0.1 and 
dashed curves for np/rieo — 0.5. We have smaller d\nT/dz 
than that for the jet, indicating a decrease in the de- 
celeration efficiency. A relatively weaker radiation drag on 
jet containing protons can also be seen in Figure 8. The 
jets with Up/rieo = 0.1 and 0.5 have much higher terminal 
Lorentz factors than the pure jet. 




Figure 8. Deceleration of jets. Deceleration in the KN regime 
and Thomson regime are indicated respectively by dashed and 
solid part of the curves. The jet is assumed to start with F = 10''' 
at 2 = 5-Rg and z = 10-Rg. In each case, the curves from left to 
right correspond to Up/n^Q = (a pure electron-positron jet), 
Wp/"eO = 0.1, and np/rieo = 0.5. The efficiency of disc emission, 
r] = 0.05, is used. Protons are assumed to be cold. 



3.3 Energy loss due to inverse Compton 
scattering 

In the calculation of the radiation force, we assume that 
the plasma in the jet is highly relativistic. In practice, the 
particle distribution is determined by the various energy loss 
processes and the injection rate of high energy particles. The 
energy loss rate due to inverse Compton scattering can be 
derived from —7 = —(3 ■ dp/dt with dp/dt given by (2.5). In 
the Thomson limit, the loss rate is obtained in the jet frame 
{Kj) as 



(l-/36Cos6ij{) -^7 



2 cRg j RO I ^ 

3Mp»?^( — ) r 7 



(3.12) 



where (...)a represents an average over Oe, and the approx- 
imation is valid only for F 3> z/Ro- In the large angle ap- 
proximation, we have —(7)0 oc 7^F^. In the small angle ap- 
proximation ^ F~^, i.e. the photon field is from a point 
source, we have —(7)0 oc 7^/F^. Prom (3.12), the time scale 
in K for energy loss for F7 < l/e^j and z > Ro is estimated 
to be 



'■ cRa 



1 

F7r? 



=0.017s 



(^) 



6 X 10-s 



F7 



10R„ 



(3.13) 



where = 

r7 > 
limit. 



6 X 10 corresponds to 30 cV photons. For 
the time scale should be calculated in the KN 
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In the KN scattering regime, their energy loss rate is 
much slower, i.e. 



the characteristic transverse size of the jet, Rj 



Ra = 1.5 X 



-{i)o 



3 Mp*? cRg , /r, ^ T^- N 



(3.14) 



Unlike (3.12) in the Thomson limit, the energy loss rate 
increases only slowly with increasing T and 7. In K, the 
characteristic time for the energy loss is 



6 X 10-5 



\lQRgJ ln(3.3r7eH) 



(3.15) 



where '^'^e^ > 1. 

Plasmas in the jet are likely to be magnetized. Then, 
particles would also radiate synchrotron emission. The en- 
ergy loss rate due to synchrotron emission is given by 



d7 



1 (TtS^7^ 

An rrieC 



(3.16) 



where B is the magnetic field in the jet frame. Thus, the 
synchrotron time scale is 



ts 



(3.17) 



The energy loss due to synchrotron emission can be compar 
rable to (3.13) and (3.15) only for strong magnetic fields or 
large 7. 

Because of the short time scales indicated by (3.13) 
and (3.15) near the black hole, for plasmas in the cell to 
be highly rclativistic one needs an efficient injection of elec- 
tron positron pairs, e.g. through cascades or acceleration. 
The relevant processes, either pair cascades or direct ac- 
celeration of electrons (positrons), must be fast enough to 
overcome energy loss due to inverse Compton scattering or 
synchrotron radiation. Although details of acceleration near 
a black hole are not well understood, the existence of a rapid 
aeecleration mechanism that is able to accelerate particles 
to ultrarelativistic energies within a short time scale cannot 
be ruled out. One possibility is acceleration by a rotation- 
induced potential drop, similar to that occuring in pulsars, 
where particles extracted from the neutron star are acceler- 
ated to ultrarelativistic energies within a much shorter time 
scale than the rotation period (e.g. Michel 1987). Particles 
injected at radial distance J?inj to the black hole are accler- 
atcd to energy 7 ~ a'^^^ AR/ Rs_n\ after travelling a distance 
Ai?, where a = eA$/meC^, A$ is the potential drop which 
can be estimated from the total power output, Lo, of the 
AGN, and AR/Rinj < 1 (e.g. Michel 1987). Then, the ac- 
celeration time tacc = AR/c is estimated to be 



CO-2/3 



10" 



''(100(22)' 



(3.18) 



where cr^/'^ ^ 10^° for Lo = lO^'^rgs'^ Since tacc < iic, 
^KN 1 electrons (positrons) can be accelerated to the energies 
in the KN regime. 

For magnetized plasmas, the condition for relativistic 
plasmas to be contained in the jet requires that the gy- 
roradius being less than the characteristic transverse size, 
Rj, of the jet. Since the gyroradius is given by pg = 
2 X 10^ cm (7/10*) (1 G/B) where B is the magnetic field, for 



lO^^m(M/lO®M0), the condition pg < Rj should be easily 
satisfied. It can also be shown that the approximation made 
earlier (Eq. 2.5), i.e. neglecting electron gyration, is valid. 



The cyclotron time scale tc ~ Pg/c 



10- 



K7/10")(1G/B), 



is much shorter than tic, tKN, ts, and even much shorter 
than the characteristic time of deceleration, which is typi- 
cally > 1 s. 

4 CONCLUSIONS AND DISCUSSION 

Radiative deceleration of relativistic jets due to inverse 
Compton scattering is discussed in both the Thomson and 
KN scattering regimes. We show that KN scattering is im- 
portant in deceleration of ultrarelativistic jets with initial 
bulk Lorentz F > 10^. For F » 1, near the black hole and in 
the jet firame, the incoming photons are beamed towards the 
hole, and so the main contribution to the drag of the jet is 
from head-on scattering by forward moving particles. These 
particles are more likely to be in the KN regime as they see 
photons of the highest energy (e' > 1). Particles scattering in 
the Thomson regime may contribute positively to the radia- 
tion force but overall KN scattering is dominant and results 
in deceleration. Thus, scattering in the KN regime should be 
included in calculations of the radiation force when F > 10^. 
Our result shows that in the KN regime, radiation drag is re- 
duced, but still severely constrains the speed of the jet bulk 
flow. Thus, Compton drag can decelerate an jet starting 
with any F > 10 in the region sufficiently close to the black 
hole down to the value F < 10. The efficiency of decelerar- 
tion is significantly reduced if the jet contains protons since 
protons have much smaller scattering cross section and the 
effect of protons is to increase the inertia of the jet. 

According to the unified scheme (e.g. Urry & Padovani 
1995), blazars are radio-loud quasars, and if this is true, 
a relativistic jet with sufficient particle kinetic luminosity 
is required to power radio lobes at a larger distance. The 
results given here show that because of radiation drag the 
ability of an electron-positron jet to power large scale radio 
lobes is severely constrained. For example, consider an jet 
with an initial luminosity NiVi ~ Ld, where Ni and F^ ^ 1 
arc the initial particle flux and Lorentz factor. As a result 
of Compton drag, the flow slows down to F/ < 10 with 
most of the energy going into radiation. Thus, even with 
pair cascades Nf 3> Ni, we may still have NfFj <^ NiVi. 
The constraint can be overcome if are re-accelerated at 
a distance further away from the central region or if there 
is outflow of accelerated protons as suggested in a two-flow 
model by Sol, Pelletier & Asseo (1989). In their model, radio 
lobes are powered by outflow of p-e plasma jets and jets 
are dominant only in the subparsec region. 

In our discussion, the external soft photons are assumed 
to come only from the disc emission which is the sum of 
blackbody emission from series of rings centred at the black 
hole, and is similar to the model discussed by Demer & 
Schlickeiser (1993). Generalization of the calculation to in- 
clude other radiation fields, such as that due to a disc torus 
(Protheroe &: Biermann 1997), and reprocessed radiation, 
should be straightforward. The effect on the jet due to scat- 
tering of photons from the broad-line region may also be 
important, and is not considered here. In calculating the av- 
erage force, we assumed an isotropic (angular) distribution 
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of electrons (or positrons) with a power law. Some mecha- 
nism, e.g. pitch angle scattering, is required for isotropiza- 
tion since inverse Compton scattering itself tends to make 
the distribution anisotropic. 
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APPENDIX A: MOMEMTUM TRANSFER DUE 
TO COMPTON SCATTERING 

Calculation of momemtum transfer (2.5) requires the full 
KN differential cross da/desdQs which was discussed in de- 
tails by Blumenthal & Gould (1970). Here, wo use the fol- 
lowing approximation (in Kj) (Reynolds 1982; Dermer & 
Schlickeiser 1993) 

« a^5[es - 7'e(l - /JcosG)] 5(0. - Qe), (Al) 
for the Thomson regime (e' < 1) and 



da 



dSsdVls 



(A2) 



for the KN regime (e' > 1), where a,^ is the Thomson cross 
section, e' is the incoming photon energy iu the electron 
rest frame, cosG = cos6cos6'e + sin 6 sin cos (0 — <j)e), Q 
is the angle between the incoming photon direction and the 
electron motion, the delta function 5{Q,s — Og) describes the 



beaming approximation. 

To evaluate (2.5), particles with 7 can be broadly sep- 
arated into two groups: one group of particles with cos 9e < 
fic) scatter in the Thomson regime with (Al) and the other 
group with cos^e > /ic scatter in the KN regime with (A2). 
This approximation is similar to that was used Dermer & 
Schlickeiser (1993) in calculating scattered photon spectrum. 
Thus, using (Al) and (A2), the radiation force (2.5) can be 
written approximately into the form 



dQed'y ~^"e(7' 



RF{R) 
27r(Ji2 -I- ^2) 



X£^r"(l -/36Cos6(^)' / dned7ne(7,ne) 



x<^ (1-/3 cos e„) 



+1 V. - " (7C0S 



r7£R 



7 (1 — /3C0S©) COS^e — COS^^ 

Efl COS 6) )//(£„ - 1; 



^(1-4) 

(A3) 



where H{e'^-\)^\ for e'^ > 1 and //(e'^ - 1) = for s'^ < 
1, ne(7, fie) is the paxticle distribution as defined earlier (Eq. 
2.7). The terms with H{1 - s'^) and H{£'^ - 1) correspond 
respectively to scattering in the Thomson and KN regimes. 

Accordingly, is written as a sum of two components, 
fz = /z,T + /z,KN where fz,T, /z.kn describe the contribu- 
tions to the force from scattering in the Thomson and KN 
regime, respectively, and they are 



dR , ^fS^K^ Snr^ii - i3b cos e^f 



2n{R^ + 

'^"dcos 9e 



ne(7,ne)(l-/3cose«) (A4) 



7 (1— /3cosQ)cos^e — COS0^ 



and 

/z,KN 



dR 



e^jF (1 — /Jfccosfe'^) 



27r(i?2 + 22) 

^dcos 



ne(7,fie)o 



3 ln(2^^4) 



(A5) 



r7£„ 



X (7 COS 0e — Efl COS ) 



where /ic, '^e(7,f2e) are defined by (2.12), (2.7). Integra- 
tion over cos^e can be easily done using approximation 



1, cosdj^ ~ —1. In this approxima- 



tion, the integrands in (A4) and (A5) can be regarded as 
independent of (j)e- Then, integrating over cos^e, we obtain 



RFe„T^ 



■ 27r(i?2 



x{7' [i(Mc - 1) - ll3 cose + 1) + i(/i^ - 1) 
-cos0^(Mc + l) + i(At?-l)}, (A6) 
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/RF 1 f 



{(l-^?)ln2- 



1 - Ate 



2£„ COS 
7 



P cos 

(l-Mc)ln(2/v^) 



(A7) 



Eq. (A6) reduces to Eq. (2.6) or (2.8) by assuming /Uc = 1, 
i.e. all particles arc in the Thomson scattering regime. Let 
It and Jkn represent 7-integration in (A6) and (A7); It and 
Ikn are calculated to be 



It = 



2e2 



7m ax 7c 



2(p + 1)4 

1 

^KN = -^(7c"''-7iir) 



3(3 -p) 



(7r"-7^r) (A8) 



1 

+- 
P 



P 

ln(4e) 



£r 



+ 2£R ln(4/e) + Aer/p 



p+l 



-^+ln(4/e) 



^O'c Tmax ^ 



(A9) 



7c ^ ln(2£R7c) - 7^f^ ln(2£H7„^J 



where 7c = rnax{l/e„ , 7^^;^}. Then, (A6) and (A7) can be 
further written into the form 



1 , - RFe^T\l-P,cos~eJ' 
^a.Coneo/ dR 27r(i?2 + ^2) 



/■Rma: 



dR 



RFIkn 



27r(i?2 + «2) 4£^ 



(AlO) 



(All) 
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